Decoupling and coherent plasma oscillations around last scattering 
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Coherent properties of the baryon-photon fluid decoupling are considered in the terms of an 
effective nonlinear Schrodinger equation for a macroscopic wave function that specifies the index of 
the coherent state. Generation of a transitional acoustic turbulence preceding formation of large- 
scale condensate in the plasma and its influence on the CMB power spectrum has been studied. A 
scaling k~ law is derived for the CMB Doppler spectrum E{k) (angle-averaged) in the wavenumber 
space, for sufficiently large wavenumber k and for the weak nonlinear and completely disordered 
initial conditions. Using the recent WMAP data it is shown that the so-called first acoustic peak 
represents (in a compensated spectral form) a pre-condensate fraction of the spectrum E{k) at a 
rather advance stage of the condensate formation process. 

PACS numbers: 98.80.Bp, 98.65.Dx, 98.70.Vc, 52.35.Ra 



I. INTRODUCTION 



The recent high resolution cosmic microwave back- 
ground (CMB) radiation measurements (see, for instance 
[l| ) provide a possibility for quantitative investigation of 
the nontrivial processes taking place at the decoupling of 
the baryon-photon plasma (fluid) in the early universe. 
Before the decoupling the CMB photons were tightly 
coupled to the baryons by photon scattering (mainly 
from the free electrons) forming so-called baryon-photon 
fluid. At a certain stage of the early universe expan- 
sion the energy of the CMB photons becomes insufficient 
to keep hydrogen ionized and a recombination process 
becomes dominating in the baryon-photon fluid. The 
recombination is a gradual process, but it serves as a 
trigger and a catalyst of the decoupling of the baryon- 
photon fluid. Before the decoupling the baryons and 
photons were tightly coupled whereas after the decou- 
pling the photons are effectively free. Physical processes 
at the decoupling itself can be described on different lev- 
els. In present paper we will concentrate on the coherent 
properties of the decoupling. It is well known that the 
main energy of the plasma fluctuations at a certain (ad- 
vance) stage of the decoupling process is concentrated in 
the coherent states (the acoustic coherent oscillations - 
stochastic standing waves). Originally the acoustic os- 
cillations become phase-coherent when the causal hori- 
zon overtakes their wavelength. Around the last scat- 
tering surface the radiation density and ionization are 
high enough that the photon drag on the baryons can be 
still significant. Therefore, the baryons are involved in 
the coherent motion by the photon drag, which differen- 
tially accelerates the ions and electrons. The resulting 
significant electron-ion drift generates large-scale coher- 
ent electric currents and fields, which rather quickly be- 
come comparable to the photon drag in their effect on 
the electrons (see, for instance, 0,[3]). While we can 
readily introduce the photon (Thomson) drag in the dy- 
namical plasma equations it is rather nontrivial task to 
describe the coherent (collective) phenomena at the de- 
coupling, which are crucial for the entire process. To 



tackle this problem we suggest to use an effective non- 
linear Schrodinger equation for macroscopic wave func- 
tion, which is a complex-valued classical field that spec- 
ifies the index of the coherent state. To the leading or- 
der each coherent state evolves along its 'classical' tra- 
jectory which is given by the corresponding nonlinear 
Schrodinger equation. The idea is based on the Ehren- 
fest's theorem stating that the the expectation values of 
displacement and momentum given by the Schrodinger 
equation obey time evolution equations which are anal- 
ogous to 'classical' ones in the the Madelung representa- 
tion. In this representation the macroscopic wave func- 
tion of the free (non recombined) electrons can be de- 
fined a,s ip = -\/'^e exp[iS'(x,i)//ic], where Ue is the elec- 
tron number density, VS'(x, t) = 777,v(x, t) with m and 
v(x, t) as the mass and velocity of the electrons, he here 
is not the Plank constant but a parameter characterizing 
coherence in the system (real quantum effects are negli- 
gible in the considered situation, cf recent papers [j| , Q 
and references therein). Then, momentum equation for 
the electron can be considered as a nonlinear Schrodinger 
equation for the macroscopic wave function ■(/'(x, t): 
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where 77(x) is a real-valued potential, C/e//(|'0P) is an ef- 
fective potential represented by a real-valued function of 
the free electrons density n^ — |V'P- Analogous approach 
at rather different circumstances was used for a mean 
field description of dense Fermi plasmas (see, for instance, 
Refs. [6], [71, 81). Unlike the quantum plasma considered 
in Refs. ^(|,[7[,[8| the original (input) coherence at the 
decoupling comes from an external source (cf Section VI). 
At this description the addition of the effective potential 
^e/zdV'P) to the equation (1) is possible due to a spe- 
cific gauge invariance of the equation. Namely, it can be 
readily shown that addition of the potential [/e//(|'0P) 
(satisfying condition lim„_,o / C^e//("-)'^"- = 0) does not 
affect the 'classic' limit (i.e. the Ehrenfest's theorem) 



md(r)/dt = (p), d{p)/dt = -{Vrj) 



(2) 



where p = mv is momentum. The effective poten- 
tial can be also directly depending on x, but the po- 
tential C/e//(x, IV'P) should be slowly varying on x (in 
the adiabatic terms of the Ehrenfest's theorem). About 
a possibility of nonlocal dependence of f/e//(|V'P) on 
IV'P see Section VI. It will be shown below that such 
type of invariance of the Schrodinger equation is cru- 
cial for the condensation process. We deliberately have 
the dissipation term not included into the equation be- 
cause we will concentrate on the large- and intermediate- 
scale processes, while the dissipation term is usually re- 
lated to the small scales (we will return to this prob- 
lem below). The first term in the right-hand side of 
Eq. (1) provides, in particular, a possibility of tunnel- 
ing at the decoupling. On the other hand, it is believed 
that the nonlinear Schrodinger equation can describe de- 
velopment of the coherent regimes @-[l2| at a certain 
stage of evolution. Theoretical considerations (see for 
instance, [9], 10]) show that this regime can appear in 
a long-wavelength region of wavenumber space after the 
breakdown of the regime of weak turbulence: the Kagan- 
Svistunov (KS-) scenario. This is supported by recent 
numerical simulations for the weakly interacting Bose gas 
[12.] and for the dense Fermi plasmas 3. It is expected 
that there exist three different stages in the time develop- 
ment for the KS-scenario: weak turbulence, transitional 
turbulence and condensate. A pre-condensate begins its 
formation in a long-wavelength region at the transitional 
turbulent regime. This transitional turbulent regime is 
the most difficult for theoretical description 12]. There 
is a problem to define a continuous vorticity filed in the 
quantum-like turbulence due to non-rotational nature of 
the velocity field defined on the wave function given by 
the Schrodinger equation. Velocity field defined on the 
wave functions is a potential field. Vorticity of such field 
vanishes everywhere in a single-connected region. It is 
believed that in the condensate state itself all rotational 
fiow is carried by quantized vortices (the circulation of 
the velocity around the core of such vortex is quantized) . 
This idea turned out to be very fruitful, and recent exper- 
iments give direct support for existence of such vortices 
in the condensate [IjI j [3 ' [3 ■ However, for the transi- 
tional turbulence apparent absence of a continuous field 
associated with vorticity is a difficult problem. Moreover, 
it is clear that the only fiux of the particles from the re- 
gion of lager energies toward the future condensate in not 
sufficient for transformation of the completely disordered 
initial weak turbulence into the condensate with its local 
superfiuid order and with the tangles of quantized vor- 
tices. A certain vorticity-like quantity have to be involved 
in this process just before appearance of the tangles of 
quantized vortices. In present paper we study a general- 
ized vorticity defined on the weighted velocity field. The 
generalized vorticity is not vanishing in the bulk of the 
fiow. Then, we show that adiabatic invariance of an en- 
strophy (mean squared generalized vorticity) controls the 
transitional turbulence just before formation of the con- 
densate with its tangles of quantized vortices. Since in 



the above described scenario the transitional turbulence 
already corresponds to a coherent regime (in this regime 
the phases of the complex amplitudes of the field tp be- 
come strongly correlated) we will speculate that an ad- 
vance stage of the decoupling can be associated with the 
transitional turbulence regime of nonlinear Schrodinger 
equation (1) . Them we will consider significance of the 
generalized vorticity for the CMB photons visibility func- 
tion, which is the main CMB characteristic of the decou- 
pling stage. This scenario provokes an additional spec- 
ulation: to associate the condensate regime itself (with 
its tangles of quantized vortices - the topological defects) 
with the universe state just after the decoupling. The 
recombination process considerably reduces the number 
of the free baryons before the last stage of the decou- 
pling but a sufficient number of the free baryons can still 
exist at this stage in order to form a condensate due to 
the fast character of the condensation process (Section 
VII and Fig. 5). The far reaching consequences of such 
additional speculation are quiet obvious, cf for instance 

[ii,[ia,[ii. 



II. 



GENERALIZED VORTICITY 



Let us recall certain well known properties of velocity 
field defined on the wave function. In order to calculate 
the expectation value of velocity one can evaluate the 
time derivative of space displacement (r) : 
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provided by the Eq. (1) one obtains 
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where we have integrated by parts with corresponding 
vanishing or periodic boundary conditions. The inte- 
grand in the right-hand-side of Eq. (4) is the rea/- valued 
distribution of the quantum-like velocity 



^ = ^1'^^'^* -^*^^^ 
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The real-valued quantum-like velocity filed itself one 
can obtain comparing the definition (v) = Jvl^Aj^dx 
with Eq. (4), that results in 
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i.e. the quantum-like velocity (6) is a potential field and 
curl V = everywhere in a single-connected region. 

Let us now define generalized vorticity on the weighted 
velocity field as 
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First term in the right-hand-side of Eq. (7) is the ordi- 
nary vorticity. When there is no the quantized vortices 
this term is zero due to Eq. (6). In this situation the 
second term in the right-had-side of Eq. (7) determines 
the generahzed vorticity. In this case lo is orthogonal to 
the "plane of motion", which is stretched over the vec- 
tors: V and Vl^/ip, at any point of the space. In this 
sense the transitional turbulence can be considered as a 
locally two-dimensional one. 

The second (linear) term r](pi)ip in the right-hand-side 
of Eq. (1) does not affect our further results and will be 
omitted for simplicity. Dynamical equation for the gen- 
eralized vorticity similarly to the dynamic equation for 
the density j-i/'P (in a dimensionless form, for simplicity) 
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does not contain the nonlinear terms explicitly 



dt 



[Vip X V(AV/)] + [Vi/;* X V(AV')] (9) 



(the right-hand-side terms in the Eqs. (8), (9) come from 
the linear term of the Eq. (1) only). This cancellation 
of the nonlinear terms in the dynamical equation for the 
density, Eq. (8), results (after integration by parts) in 
the conservation law for the total number of the particles: 
N = J \ip\'^dx = const. The cancellation of the nonlinear 
terms in the dynamical equation for vorticity, Eq. (9), 
results in a less strong result. Namely, the enstrophy 
(mean squared generalized vorticity) turned out to be an 
adiabatic invariant for the transitional turbulence. 



III. ADIABATIC INVARIANCE AND SCALING 
TRANSITIONAL REGIME 

According to the theoretical predictions [1] , [l3| if one 
starts from a self-similar solution of the equation (1) for 
weak- nonlinear conditions, then the so-called coherent 
regime 0,[l3|,[ll|, will be developed at a certain stage 
of evolution. The first stage of evolution leads to an 
explosive increase of occupation numbers in the long- 
wavelength region of wavenumber space where the order- 
ing process takes place. From the beginning of the coher- 
ent stage of the evolution we will call the long- wavelength 
region of wavenumber space as pre-condensate fraction, 
while the rest of the wavenumber space we will call as 
above-condensate fraction. 

At a certain time, close to the blow-up time of the self- 
similar solution, the coherent regime sets in. After this 
time the system has a certain transitional turbulent pe- 
riod. At the end of this transitional period appearance of 
a well-defined tangle of quantized vortices indicates the 
final (condensate) stage of the evolution. Therefore, gen- 
eralized vorticity exchange between the pre-condensate 
and above-condensate fractions in the transitional period 
seems to be crucial for the process of the condensate for- 
mation in the KS-scenario. 



In the case of space isotropy it is convenient to deal 
with an angle-averaged occupation numbers spectrum Nk 
in variable k = |k| in the wavenumber space 



N == / Nkdk 



(10) 



The total number of particles associated with the pre- 
condensate fraction is N^'^ = Ik'<k ^'k^k', where fcc 
is the wavenumber scale separating (approximately) the 
pre- and above-condensate fractions (for its determina- 
tion from numerical simulations see next section). Then 
the total number of particles associated with the above- 
condensate fraction is: N"-'' — N ~ N^"^. From the dy- 
namical conservation law for the total number of parti- 
cles, N — const, we obtain 



dNP'= 
dt 



dN" 
dt 



(11) 



Then we can define exchange rate of the particles number, 
e„ and as 



(12) 
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Let us also consider angle-averaged enstrophy spec- 
trum life 



(^') = y"n, 



dk 
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Then, the enstrophy associated with the pre-condensate 
fraction is ^^'^ = Jk'<k H^^'' whereas enstrophy as- 
sociated with the above-condensate fraction is: 0.°''^ = 
{lo'^) — VLP"^. The intensive particle flux to the long- 
wavelength region of wavenumber space (where the pre- 
condensate is formed) involves the enstrophy one. Unlike 
the total number of the particles, the average enstrophy 
{uj^) is not an exact invariant of the motion. Though, 
characteristic time scale of the enstrophy exchange be- 
tween the pre- and above- condensate fractions is ex- 
pected to be much smaller than the characteristic time 
scale of the (w^) evolution Q-IHI. Therefore, the (w^) 
can still be considered (approximately) as an 'adiabatic 
integral' for the exchange process. As it is used for the 
adiabatic processes this statement can be formalized as 
following: 
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\dn''''/dt\ d\{uj'^)/dty \dnp''/dt\ \d{uj^)/dt\ 

(14) 
(actually, only one of the (14) inequalities is sufficient 
for further consideration). For the nonlinear system the 
principal problem is invariance of the characteristic time- 
scales in (14) to the rescaling ip — > Xip (otherwise the 
inequalities (14) become meaningless). It can be readily 
shown that the cancellation of the nonlinear terms in the 
dynamical equations (8), (9) provides such invariance for 
the enstrophy based time-scales in Eq. (14). 



Since for the transitional turbulence fi"'^, ri^"^, and {uP') 
are still of the same order, then we obtain from Eq. (14) 
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Then from n"-" = {uj'^) - nP" and Eq. (15) we obtain 
dnP'= dVL"" 



dt 



dt 



(16) 



and we can define exchange rate of the enstrophy, e^} in 
full analogy with Eq. (12) 
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i.e. the adiabatic invariance of the enstrophy substitutes 
its exact invariance in the case of the transitional turbu- 
lence. 

For the weak nonlinearity one can expect that the 
angle-averaged spectrum Nk for sufficiently large k is pro- 
portional to e„ (cf [24|,l23|). This can be also valid for 
the transitional turbulence for sufhciently large k. Then 
with the additional dimensional parameter e^ (which has 
dimension T~^) one obtains scaling law for sufficiently 
large k from dimensional considerations 
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(18) 



The enstrophy exchange controls the particle exchange 
in this case. 

External dissipation or forcing can make the estimates 
(14) invalid and, hence, destroy this scenario (cf [24|.[25|). 
Though, in the case of a sufhciently weak and linear ex- 
ternal dissipation (forcing) the above consideration can 
be still valid. In this case the exact invariance of the total 
number of particles N can be replaced by its adiabatic 
invariance (similarly to the enstrophy) and the inequali- 
ties similar to the Eq.(14) can be made. Due to linearity 
of the dissipation (forcing) these inequalities will be still 
invariant to the rescaling: ijj — > Xip. 

There exists another way to tackle the problem. The 
wave function can be considered as a vector on the com- 
plex plane. The density: \tp\'^, gives magnitude of the 
vector, whereas the phase gives its direction. Conditional 
functional average over the directions of the vector can 
reduce the dynamics of the |V'P to a passive scalar one 
(cf [30|). Then, spectrum like Eq. (18) can be obtained 
in analogy with classic fluid turbulence controlled by the 
enstrophy adiabatic invariance (cf [3l|). 



IV. COMPARISON WITH NUMERICAL 
SIMULATIONS 

In paper |12l | a large scale t/iree-dimensional numer- 
ical simulations of the Gross-Pitaevskii equation: with 
r7(x) = and C/e//(|V'P) oc IV'Pj was performed in order 
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FIG. 1: The cumulative number of particles iV^ — Efe/<fc A^^/ 
against Infc. The data are taken from [12t | for t = 600 (the 
beginning of the transitional stage) . The straight line is drawn 
in order to indicate agreement with Eq. (19). 



to reveal all three stages of the evolution from weak tur- 
bulence to superfluid turbulence (the Bose-Einstein con- 
densate) with a tangle of quantized vortices. While the 
observed evolution for t < 600 exhibits the well defined 
self-similar weak turbulence, the period 600 < t < 1000 
was identified in |12l | with the transitional turbulence. In 
Figs. 1 and 2 one can see the cumulative number of par- 
ticles Nk = Sfe'<fe ./V^, which shows how many particles 
have momenta not exceeding k at t = 600 (the beginning 
of the transitional stage. Fig. 1) and at t — 1000 (the 
end of the transitional stage. Fig. 2). 

Using scaling (18) one can estimate the number for 
sufhciently large k as 

Nk = Sfe'<fe Nk' - Sfe^<fe'<fc [k')-^ ~ \n{k/kc) (19) 

where kc is the wavenumber scale separating (approxi- 
mately) the pre- and above-condensate fractions. We use 
the semi-log scales and the straight lines in the Figs. 1,2 
in order to indicate agreement with the Eq. (19). After 
the formation of the quasi-condensate at t > 1000 (third 
stage), the distribution of particles acquires a bimodal 
shape [12l | . The value kc can be defined at crossing point 
of the straight line in Figs. 1,2 (indicating Eq. (19)) and 
the horizontal axis. 



V. VISIBILITY FUNCTION AND CMB 

So-called visibility function is used to describe 'optical' 
properties of the plasma at decoupling stage. The visi- 
bility function g{t) is defined so that the probability that 
a CMB photon last scattered between time t and t + dt is 
given by P{t)dt. A measure of the width of the visibility 
function around its (rather sharp) maximum can be used 
for quantitative characterization of the decoupling stage. 
This stage is also called as last scattering shell (or for 
the CMB measurement purposes as last scattering sur- 
face). The last scattering probability turns out to be a 
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FIG. 2: As in Fig. 1 but for t = 1000 (the end of the transi- 
tional stage). 
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FIG. 3: The cumulative compensated spectrum of the CMB 
temperature fluctuations calculated using the WMAP 3-year 
data [3|. The dash straight line in this figure indicates agree- 
ment with Eq. (27). The cumulative spectrum is given in ar- 
bitrary units. The insert shows the WMAP three-year binned 
spectrum l{l + 1)/2tvCi against I [l| (see Eq. (32)). 



narrow peak around a decoupling redshift. The visibility 
function is defined as 



ne{t')aTa{t')dt' \ (20) 



g(t) — ne<JTa{t) exp 



where a{t) is the expansion factor normalized to unity 
today, Tie is the electron density, and ax is the Thom- 
son cross section. The motion of the scatters imprints 
a temperature fluctuation, 5T, on the CMB through the 
Doppler effect 



5T{n) 



g{L)n ■ v(x) dL 



(21) 



where n is the direction (the unit vector) on the sky and 
V is the velocity field of the electrons evaluated along the 
line of sight, x = Ln. Using Eq. (20) let us introduce a 
modified visibility function: g' ~ gjn^. Then, 
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where we used Ue = |V'P and j = IV'pv is the weighted 
velocity filed. Unlike the potential velocity field v (Eq. 
(6)) the weighted velocity j has a rotational component 
in addition to the potential one (V x j ^0, Eq. (7)). It is 
significant in present content due to the well known 'geo- 
metrical' cancellation of the contribution of the potential 
component into the integral (22) (see [34|,[33]). The well 
known Helmholtz theorem decomposes any sufficiently 
smooth, decaying vector field into potential (curl-free) 
P and solenoidal (divergence-free) S component vector 
fields: j = P + S, where 



S(r) = -^V X 
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so that 
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Since V x '^ — lo\'4>\'^ , then taking time derivative of the 
both sides Eq. (23) and using Eq. (9) one can imme- 
diately conclude that dynamical equation for S does not 
contain the nonlinear terms explicitly (similarly to the 
equations (8) and (9)). Therefore, one can define an ex- 
change rate of S^: Eg, in full analogy with Eq. (17). 
If the weak nonlinear approximation is valid for the de- 
coupling stage, then, the angle-averaged spectrum of the 
S-field fluctuations (S^) — J ^kdk: 
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can be analogously obtained for sufficiently large k. That 
means that the scaling k~^ can be also observed for corre- 
sponding angle-averaged CMB temperature fluctuations 
spectrum E{k) 
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(26) 



for sufficiently large k. It follows from Eq. (26) that the 
cumulative compensated spectrum 



k'E{k') 



k'E{k')dk' ~A + Bk 



k'<k k'<ka 

(27) 
where A and B are certain constants. The cumulative 
spectrum calculated using the WMAP 3-year combined 
data (version 2.0, March 2006) is shown in figure 3. The 
dash straight line in this figure indicates agreement with 
Eq. (27) (the wavenumber k is normalized by the co- 
moving angular-diameter distance to the last scattering 
surface: R). 



The resonant acoustic waves (see Introduction) present 
a finite-size effect for this scaling (actually, the so-called 
first acoustic peak represents the pre-condensate, see 
next Section). Therefore, these waves can be taken into 
account in a standard for the finite-size effects way [36[: 



E{k) ^ f{k/K)k- 



(28) 



where the function f{k/kc) of the dimensionless 
wavenumber k/kc corresponds to the finite-size effect. 
Then the properly compensated spectrum 



kE{k) ^ f{k/k. 



(29) 



represents the acoustic-resonant effect after extraction of 
the scaling component (26). 

In the CMB literature it is used to expand the CMB 
temperature fluctuations in spherical harmonics 
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FIG. 4: Angular-average energy spectrum for a two- 
dimensional dense plasma. The data are taken from a numer- 
ical simulation presented in [g| for the Eq. (33) with weak 
nonlinear conditions. 



For isotropic situation the angular power spectrum of the 
fluctuations, C;, is defined as {o-imO-'i' m') ~ CiSu'Smm' . 
Then, for sufficiently large I 
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(/ + l)Ci dl (31) 



(see for more details [355). It follows from Eqs. (26) and 
(31) that for sufficiently large /: E{k) (x {I + l)/27r Q, 
where the comoving wavenumber k ~ 1/ R. Hence, for 
sufficiently large I the finite-size (i.e. acoustic-resonant) 
corrected scaling E[k) ^ f{k/kc)k~^ can be written in 
the compensated form kE{k) ^ f(k/kc): 



l{l + l)Ci^f{l/Q 



(32) 



where f{l/lc) represents the acoustic-resonant peaks. It 
should be noted that in the observed angular power spec- 
trum Ci and even in the (l + l)Ci (which represents spec- 
trum E{k) for sufficiently large I ~ kR) one cannot see 
a hint on the first acoustic-resonant peak (see next Sec- 
tion). The only properly compensated spectrum kE{k) 
(Eq. (29)) or {{l + l)Ci (Eq. (32)) gives clear exhibition 
of the acoustic-resonant peaks (see the insert in Fig. 3). 
It is interesting to note that for small I [i < I < 30) 
the scale- invariant potential perturbations generate CMB 
temperature fluctuations, which also have a nearly con- 
stant l{l + l)Ci (the so-called Sachs- Wolfe effect (Hi). 
This makes the compensated spectrum l{l ^ l)Ci an ad- 
equate tool for a wide range of scales. 
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where the potential $(x, |?/jp) is a functional of the den- 
sity |^(x)p, e is the electron charge, e the dielectric con- 
stant, and no is the equilibrium electron number density. 
Splitting the potential $ = $i -I- <I>2, where 
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(34) 



and using the integral representation of the potential <I>i 



$r(x,|^n^-- 
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(35) 



it can be readily shown that the nonlinear Schrodinger- 
Poisson system possesses the same type of mesoscopic 
gauge invariance as the Schrodinger equation (1) dis- 
cussed in the Section I. I.e. similarly to the nonlin- 
ear potential C/e//(|'0P) the nonlinear nonlocal potential 
<i>i(x, \ip\^) does not appear explicitly in the 'classic' fimit 
Eq. (2) (i.e. in the Ehrenfest theorem). Indeed, 



(V$i) oc 



1 



i^(x)n^(x')pv 



|V^(x)nV;(x')P f^^dxdx' 



dxdx.' 



(36) 



VI. PRE-CONDENSATE 

For the quantum Fermi plasmas an additional (nonlo- 
cal) nonlinearity was also considered in the form of a po- 
tential given by an additional Poisson equation [y],[7|,[3|: 



dTp 



— V2 + r;(x) + <i>(x,|V'P)-C/e//(|^P) 



^, 



In recent paper [g] results of a two-dimensional numer- 
ical simulation with a dense Fermi plasma have been re- 
ported. The authors of [g| have used the equations (33) 
with Ueff{\ip\'^) (X IV'P- An energy spectrum obtained 
in this numerical simulation under weak nonlinear condi- 
tions is shown in Fig. 4 in the log-log scales. The authors 
of the Ref. Q related the strong hump in the infrared 
part of the spectrum to a condensation process, which 
was also directly observed in this simulation. Although 
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FIG. 5: Angular-average spectrum E{k) (Eq. (37)) in the 
wavenumber space for the binned WMAP data il] for I = 
kR > 30 (cf. the insert in Fig. 3). 



the causes for condensation and the coherent plasma os- 
cillations in the dense Fermi plasma considered in Q are 
rather different from those at the decoupling, the equa- 
tions and the main properties of the condensation process 
are similar. Therefore, it is interesting to compare the 
results of the simulation with the CMB data. To make 
this comparison possible one needs in the spectrum E{k) 
for the temperature fluctuations. As it is shown in the 
Section V for sufficiently large I: 



E{k)o,ikR+l)CkR 



(37) 



where in the {I + 1)C/ the variable I is taken as / ~ kR. 
Considering Z = 30 as a sufficiently large value for this 
purpose we show in Fig. 5 spectrum E{k) (in the log- 
log scales) calculated using the binned WMAP data [l| 
(these data are shown in another presentation in the in- 
sert in Fig. 3) . The arrow in Fig. 5 indicates the absolute 
maximum of the graph shown in the insert in Fig. 3 (the 
top of the first acoustic peak) . Similarity between figures 
4 and 5 is quiet obvious. It seems plausible to interpret 
the scale kcR (Fig. 5) as the scale separating between 
the pre- and above-condensate fractions at the decou- 
pling (see Sections III, IV). This is also the input-scale of 
the external coherence at the decoupling (see Section I). 
The flatness of the pre-condensate fraction of the spec- 
trum E{k) (Fig. 5) indicates a rather advance stage of 



the condensation at the end of the decoupling process (a 
flnal part of the transition to the condensate state) . 



VII. DISCUSSION 

Different sources of turbulence in the early universe 
in the times preceding the decoupling were discussed in 
the literature, from the very early times [331 up to the 
recombination time (38[-[49|. In present paper we have 
discussed a possibility that the decoupling process itself 
can generate a turbulence related to a condensation pro- 
cess in the primordial plasmas. Just this 'later' turbu- 
lence has a preferable chance to imprint itself both on the 
CMB and on the further dynamics of the luminous mat- 
ter. Especially significant is the question: whether the 
condensate itself with its tangles of quantized vortices 
appeared at the end of decoupling or the decoupling was 
finished at the transitional stage of the condensation pro- 
cess (see the end of the Section I) . The apparent flatness 
of the pre-condensate fraction of spectrum E{k) in Fig. 
5 (in the range: kcR > kR > 30) indicates a rather ad- 
vance stage of the condensation process in this case. One 
can even speculate that the end of the transition from 
the pre-condensate to the condensate is the effective end 
of the decoupling (and just this is imprinted in the CMB 
spectrum shown in Fig. 5). But the question is too se- 
rious to make definite conclusion from this observation 
only. 

On the other hand, the recombination and dissipa- 
tion processes could effect the decoupling turbulence. 
However, the coherent nature of this comparatively 
large-scale turbulence makes it stable to these effects. 
First of all, the recombination obviously is much slower 
process then the particles (enstrophy) exchange between 
the pre- and above-condensate fractions. Therefore, 
the total number of the particles N can be considered 
as an adiabatic invariant at the decoupling and the 
consideration of the section III can be preserved when 
one takes into account the recombination. The linear 
dissipation can be treated by the same way (considering 
N as an adiabatic invariant). 
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